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We present a general theory for using an optically pumped diamond nitrogen-vacancy center as 
a tunable, non-equilibrium bath to control a variety of nuclear spin dynamics (such as dephasing, 
relaxation, squeezing, polarization, etc.) and the nuclear spin noise. It opens a new avenue towards 
engineering the dissipative and collective nuclear spin evolution and solves an open problem brought 
up by the 13 C nuclear spin noise suppression experiment [E. Togan et al, Nature 478, 497 (2011)]. 
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Introduction. Diamond nitrogen-vacancy (NV) center 
11] is a leading platform for quantum computation and 
nanoscale sensing HI- The NV spin and a few sur¬ 
rounding nuclear spins form a hybrid quantum regis¬ 
ter M- Its coherence time is ultimately limited by 
the noise from environmental nuclei. This motivates 
widespread interest in using the NV spin to control the 
qubit and environmental nuclei through their hyperfine 
interaction (HFI). In addition to the remarkable success 
in manipulating fill . ITit a few qubit nuclei, there is in¬ 
creasing interest in controlling the nuclear spin dissipa¬ 
tion, e.g., dephasing and relaxation of individual qubit 
nuclei 13, [14] and dynamic polarization of many envi¬ 
ronmental nuclei [15 h 18I ]. Intensive experimental efforts 
have led to dramatic enhancement of the NMR signal 
0 for applications in chemistry and biomedicine and 
the first demonstration of coherence protection by sup¬ 
pressing the nuclear spin noise 0 - This is an important 
step towards engineering the nuclear spin evolution for 
coherence protection, nanoscale sensing [31.20—22]], and 
long-time storage of quantum information |23l . 1241 ]. 

This prospect, however, could be hindered by our lim¬ 
ited understanding of the dissipative nuclear spin dynam¬ 
ics. At present, theoretical treatments are limited to p he- 
nomenological or semiclassical modelling 13, MMM 

t or numerical simulation neglecting the NV coherence 
0- The former provides an intuitive picture, but is 
qualitative. The latter is more accurate, but is limited 
to a small number of nuclei and may miss important ef¬ 
fects due to the NV coherence. Crucially, it is not clear 
how to efficiently and quantitatively control the nuclear 
spin dissipation and especially the nuclear spin noise, e.g., 
the physical mechanism leading to the most impressive 
observation of Ref. [0 . the unconditional suppression 
of the 13 C nuclear spin noise without appreciable polar¬ 
ization, remains unclear. Subsequent noise suppression 
experiments [2(1, 27] are based on the simple but chal¬ 


lenging approach of completely polarizing all the nuclei 
or conditioned on measurement-based postselection 0 . 
Despite the recent experimental progress in controlling 
the nuclear spin polarization in certain setups 0 .30 , a 
general guidance for the efficient, unconditional control 


of the nuclear spin dynamics and noise is still lacking. 

In this letter, we present a quantum theory for using 
the NV center to engineer various nuclear spin dynamics 
and noise. The essential idea is to introduce tunable dis¬ 
sipation into the NV center by optical pumping, so the 
NV becomes a tunable , dissipative bath for the nuclei. 
When the NV dissipation is much faster than the NV- 
induced nuclei dissipation (i.e., the bath being Marko¬ 
vian), we derive a generalized Lindblad master equation 
for the many-nuclei density matrix, with analytical ex¬ 
pressions for the nuclei transition/dephasing rates. They 
not only allow easy calculation of various nuclear spin 
dynamics incorporating the NV coherence, but also al¬ 
low engineering these dynamics (dephasing, relaxation, 
squeezing, polarization, etc.) and the nuclear spin noise 
by controlling the NV. This is illustrated by (i) control of 
nuclear spin relaxation and dephasing, (ii) nuclear spin 
squeezing, and (iii) suppression of the noise from many 
13 C nuclei. Case (i) provides a microscopic basis for the 
phenomenological spin-fluctuator model 0] and experi¬ 
mental observations UM , and a simple method to sup¬ 
press the nuclear spin dephasing or relaxation. Case (iii) 
provides a general and efficient way to suppress or am¬ 
plify the nuclear spin noise and explains the observed 13 C 
nuclear spin noise suppression [19( as a special case. 

General theory.-'We consider many nuclei {Ifc} (de¬ 
scribed by the Hamiltonian Hn) coupled to an optically 
pumped NV center. We always work in a suitable NV ro¬ 
tating frame and nuclei interaction picture, so the Hamil¬ 
tonian consists of the time-independent NV part H e , the 
longitudinal HFI K that commutes with Hn, and the 
transverse HFI V ( t ) that flips the nuclei: 

p(t) = -i[H e + I< + V(t), p(t)] + C e p(t), (1) 


where C e p = 7/*^ , [I/)(*|]p is the NV dissipation in 

the Lindblad form T>[L}p = LpL t — {L^L, p}/2. Here we 
focus on NV-induced nuclei dynamics. The direct nuclei 
interactions and intrinsic nuclei damping can be easily 
included at the end of the derivation. 

To derive a closed description for the many-nuclei state 
p(t) = Tr e p(f), we use the adiabatic approximation [32- 
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34) to eliminate the fast electron motion. We define the 


many-nuclei basis |m) as the common eigenstates of Hjy 
and K with Jv|m) = A' m |m), where I\ m is an electron 
operator. The block /3 m , n = (m|p|n) obeys 

Pm,n = Pm.nPm.n l{pm,nj n }/2 ) P\ I*- 1 }? 

where 5I\ m n = A m A n and A m n (*) = z[JT e T(A m T 
A' n )/2,*] + £ e (*). Tracing over the electron yields 

p m.n = i Tr e {/5 m , n , <5A m-n }/2 i Tr e (m| [V”, /)] |n) 

for p m ,n = (m|p|n). The above equations contain three 
dissipation time scales: NV dissipation (time scale T e ) 
driven by £ mn , nuclei dephasing (time scale T 2 ) by 
5K m rl fluctuation, and nuclei relaxation (time scale T\) 
by V ( t ) fluctuation. Nuclei dissipation much slower than 
T e can be adiabatically singled out. For specificity, we 
consider T e -C T \, T 2 and single out all the dynamics of 
p(t) on the coarse grained time scale At T e . 

By treating C m n exactly and 5K m ni V(t) perturba- 
tively, application of the adiabatic approximation 


gives the nuclear spin dynamics order by order p = 
(p)i + (p ) 2 + • • •. The first-order dynamics describes nu¬ 
clear spin precession in the electron Knight fields, 


(P) 1 = ~*Tr e [K + V(t),p 0 (t)], 


( 2 ) 


equivalent to a renormalization of the nuclei Hamiltonian 
H n , where p 0 {t) = J2 m ,n |m)(n|p m , n (t).P m , n and P m ,„ 
is the normalized electron steady state: £ m , n -P m , n = 0. 
For (p|V(t)|m) = Fe~ lbJt , we obtain 


(Pm,m)2 — pPp,p ^pHnPm,m)! 

e iut Tr P A t e £p ' m * FP m n 


P 


Wp 4 _ m = 2 Re 


dt. 


( 3 ) 

( 4 ) 


The off-diagonal coherence p m ,n (m / n) obeys 

^ _ (-rip , Sp(hbp < -n|m + 1 (Tp<-m| n) 

(Pm,nj2 — ( Y + 


Pm.n 1 

( 5 ) 


= Re 


Tr e 5K n 


at SK m ^P mtn dt, ( 6 ) 


where SK mn = d /\ ln n — Tr e SK , 
for Wp 4 _ m | n is involved 


m.nAn.n- The expression 


but it reduces to W, 


p«— m 

upon neglecting the difference between K m and K n . The 
key quantities of our theory, the transition rate Wp^ m 
[Eq. ©] and pure depliasing rate T^ n [Eq. ©] are 
obtained by calculating the inverse (£ p . m + iw ) _1 and 
^-ain — hm^_s.o(A m .n + iv)~ l . Here we notice that if 
we focus on the dynamics of {p m . m} on the time scale 
At T e ,T 2 , then we can treat SK mn exactly and still 
derive Eqs. © and with £ p . m in Eq. © replaced 

with £d> ot m = Ap, m (*) - *{•, <5A'p,m}/2. 


Equations m describe a variety of dissipative and 
collective nuclear spin dynamics governed by the NV. 
Before engineering them, we present a perturbative ex¬ 
pression for Wp<_ m to exemplify the previously ne¬ 
glected effect of NV coherence. For simplicity we 
set AT = 0, so £ m , n = A and P m , n = P are in¬ 
dependent of the nuclei state. For (p|V(t)|m) = 

Vfie~ lult \f){i\ (/ ^ i), Wp<_ m consists of the golden rule 
part 2'ir\Vfi\' 1 {i\P\i)5^ Tt+Vi P 2 ' , (Y{,ezf.i) and the coherent 
part W^ m = 2 \Vf t i\ 2 ImJ2j&{ i \P\j){j\He\i)/{zf, i Zfj), 
where Tj = YlflfU ^\x) = {^/tt)/(x 2 + y 2 ), and 

z k ,j = (k\H e \k) — (j\H e \j) —uj — i(Tk + Tj — 2"/k,k5k,j)/‘2' is 
the complex energy mismatch. Under optical/nricrowave 
driving, NV coherence and Wp^ m could be important, 
e.g., it dominates 13 C nuclei flip by the NV ground state 
in the first nuclear spin noise suppression experiment fljj | 
(to be discussed shortly). For more general (p|V(i)|m), 
an explicit expression for W pi - m is given in [35[. It can be 
easily used to calculate the nuclear spin transition rate in 
a given experimental setup. Now we illustrate controlling 
the nuclear spin evolution by manipulating the NV. 

Nuclear spin dephasing and relaxation. -We consider a 
single 13 C or 14 N nucleus upon exciting a cyclic NV op¬ 
tical transition | g) -H- |e) with detuning A, e.g., | g) = |0) 
and |e) = |&,) in the widely used setup for single-shot 
readout [2a, [37). Dropping the flip between different 
NV states, the HFI takes the general form F • I, with 
F = \g)(g\a g + |e)(e|a e the NV Knight field. The NV dis¬ 
sipation includes the radiative decay (rate 71 ) from |e) to 
| g) and the pure dephasing (rate 7 ^) of |e). Equation © 
gives (p) 1 = — i[F ■ I,p], where the averaged Knight field 
F = TrFP defines the nuclear spin z-axis e z = F/|F|, P 
is the steady NV state with (e|P|e) = W/( 7 i + 2 W) and 
W = 27r(Di?/2) 2 (i^ 7l+7 'A/ 2 )(A). For |a g |, |a e | <C 71 , Eqs. 
©-© give the second-order dynamics in the Lindblad 
form: (p) 2 = 2Y v T>[i z ]p+ Y 1 (V[I + \ +2?[/_])p, where 

© = co(e|P|e)( 7 i + 2 W)~\a; - a*) 2 , (7) 

Fr = (co/2)(e|P|e)( 7 i + 2W)~ 1 \( a s - a e ) ± | 2 , ( 8 ) 

and Co is a dimensionless 0(1) quantity. The average 
nuclear spin obeys d t (I z ) = — (A)/Ti and d t {I+) = 
-(I+)/T 2 with T x = 1/(2F!) and T 2 = \/{T v + r,). 
Equation © [©] shows that the nuclear spin pure de- 
pliasing (relaxation) is controlled by the fluctuation of 
the longitudinal (transverse) Knight field. This was first 
pointed out in Ref. [13) , where a phenomenological spin- 
fluctuator model was proposed for numerical simulation. 
Our analytical results Eqs. © and © not only provide 
a microscopic basis for the previous model [13| and ex¬ 
perimental observations [2j,[3l|] (e.g., it clearly shows mo¬ 
tional narrowing r v ,Fi oc 1 /IV under saturated pump¬ 
ing W ^ 71 ), but also demonstrate the possibility (l3[] to 
control Yy and Ti by a magnetic field B: since the nuclear 
Zeeman term — 7 jvB ■ I renormalizes a g / e to a g / e — 7 jvB, 
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we can always tune the nuclear spin quantization axis 
e, oc F = (g\P\g)a g + (e|P|e)a e - 7 jvB to e z 1 a g - a e 
(e 2 || a g — a e ) such that T^, = 0 (Tj = 0). Interest¬ 
ingly, the sum rule T^, + 2T\ oc \a g — a e | 2 suggests that 
reducing T v (ri) inevitably increases Id (T v ) and it is 
impossible to suppress r v and Ti simultaneously, unless 
the NV states are tuned such that a g = a e . For more 
NV levels, analytical results are no longer available, but 
our general theory is still applicable. 


Nuclear spin squeezing.- Here we explore 13 C nuclear 
spin squeezing [381 ] by engineering the first-order evolu¬ 
tion Eq. ©• With a magnetic field Be z to quantize 
all 13 C nuclei along e, || N-V axis and a microwave 
to couple the NV ground states |0) and | + 1} with 
detuning A, the rotating frame Hamiltonian consists 
of H e = (H fl /2)(| + 1)(01 + h.c .) + A| + 1}(+1| and 
K = h\ + 1)(+1|, where h = £„<+l|S| + 1) • A n ■ I n « 
J2 n a n I ntZ comes from the dipolar HFI with 13 C nuclei 
and a n = (+1|S| + 1) -A n -e z . To introduce tunable dissi¬ 
pation, we consider weak optical pumping of | + 1 ) (with 
rate R) to the orbital excited state |e, +1), which decays 
back to | + 1) (with rate 7 ), or to a singlet | S) (with rate 
7 i c ) and then to |0) (with rate 7 S ). This creates a unidi¬ 
rectional transition from | +1) to |0). The transition rate 
71 is tunable from Ry lc /(7 + R) (small R) to 7 S (large 
R). We define the nuclear spin basis (|m)} as eigen¬ 
states of h with eigenvalues {h m }- For 71 \h m — h n |, 
Eq. © gives (p m ,„)i ~ -i(m\[hPu(h),p]\n), where 
P\i(h) = W{h)/{^i + 2W(h)) is the ^-dependent pop¬ 
ulation of | + 1) and W(h) = 2'k{VLr/2) 2 8^ 1 ^' 2 \N. + h). 
According to Ref. [38j, for polarized nuclei initially along 
e x axis (prepared by rotating e~-polarized nuclei by a r.f. 
pulse), the evolution under the nonlinear Hamiltonian 
hPu(h) ~ Pn(0)/i + P 1 ' 1 (0)ft 2 could lead to nuclear spin 
squeezing, even for non-uniform coupling {a n }. Taking 
a n = a for an estimate, the characteristic squeezing time 
38, 39] for N nuclei is = [P 1 ' 1 (0)A r a 2 ] -1 . The maxi¬ 
mal iV-nuclei collective dephasing rate ~ N 2 T V with T^ 
obtained analogous to Eq. 0 - For A = = 2yi, we 

have T^ V) 4 Ar) ~ IV/100, suggesting significant squeezing 
for N -C 100 nuclei without appreciable dephasing. 


Controlling nuclear spin noise.- Here we consider con¬ 
tinuous pumping the NV to control the noise from 
many 13 C nuclei coupled to the NV via dipolar HFI 
Y2n=i S-A n -i„. To focus on noise control, we neglect the 
squeezing effect, so Eq. © gives (p)i = b n -I„,p]. 

Here b„ = (TrSP) • A n — y^B defines a local coor¬ 
dinate (e n} x,e ntY ,e n ,z) for I„, where e n , z oc b„ and 
P is the steady NV state in the absence of 13 C nuclei. 
We decompose the dipolar HFI into k = S • h ~ S z h z 
(e z along N-V axis) and the remaining part V, where 
h = ^ n (A„ • e n,z)In,z- The I\ term not only allows 
the thermal fluctuation of h z to rapidly decohere the NV 
spin, but also allows the NV spin magnetometer to mon¬ 



FIG. 1. (color online) (a) NV center under CPT at low tem¬ 
perature 0 . Solid (Dashed) arrows denote laser excitation 
(Lindblad damping). Theoretical results (lines) vs. experi¬ 
mental data [0 (symbols) for (b) 14 N population on |mo = 0} 
and (c) NV fluorescence for different preparation of 13 C states. 
Relevant parameters: A e = 40 MHz 40|, A g = 2.2 MHz [ill ]. 
7 si = 7 = 1/(12 ns), 7 S 2 = 7 / 120 , 7 ce = 7 / 8 OO (obtained 
by fitting the fluorescence data [Tg|). -y s = 7/33 jf|, photon 
collection efficiency e = 5x 10” [iJJ], 7 c = 2.5 x 10 ” 2 s” 1 , 
Pcond = 0.288 ms, ftA = 2 MHz, readout Rabi frequency 
STa =3.2 (black), 10 (orange), and 8 MHz (blue). 


itor and control the slow h z fluctuation by engineering 
the following feedback loop [0. Here for simplicity we 
take A„ = A, so that h z = aJ2 n k,z < h max = aN/2 
and a = e z ■ A n ■ e n ,z, although the physical mecha¬ 
nism is general: (i) Through k = S z h z , the NV spin S z 
monitors the fluctuation of h z and records its instanta¬ 
neous value h into the NV steady state P(h), as deter¬ 
mined by — i[H e + S z h,P] + C e P = 0; (ii) Through V. 
the NV state P(h) flips each 13 C with rates 
and Wi^{h) [given by Eq. © with P m , m replaced 
by P{h) ] and drives h z from h to the steady-state value 
H(/i) = NaJ2 n {In,z), where (J n ,z) = {l/2)\W^k h ) - 
Wi^{h)\/\W^±{h) + W\,^{h)\. For example, to lock 
h z to a pre-defined value /i pre , we can design the feed¬ 
back loop such that HI(ft. pre ) = /i pre and the derivative 
Et'(/i P re) < 0, i.e., the NV flips the nuclei to decrease 
(increase) h z upon detecting h z > h pie (h z < h vie ). To 
describe this feedback loop, we quantify the h z noise by 
the width of the h z distribution p(h) = Tr5(h z — h)p in 
the A-nuclei steady state p and apply the nuclear spin 
feedback theory [0 to Eq. © and obtain 


p(h) oc 


e -(h-h,f/(2a 2 ) 

[W^(h) + W^ t (h)}[ 1 - hW(h)/hl ax ] ’ 


(9) 


with /i* = H(/i*) the most probable value of h z , cr 2 /< 7 2 h = 
(1 — ft 2 /^max)[l — and a t h = ay/N /2 the ther¬ 

mal fluctuation of h z . Equation © summarizes three 
ways to suppress the h z noise, (i) Narrow e~( h ~ h *) /( 2ct 1 
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by high nuclear polarization, e.g., ft*/ft max ~ ±90% re¬ 
duces a by a factor of 2. (ii) Narrow e -0- ft *) A 2a ' 
by strong negative feedback H'(ft*) <C — 1. This pro¬ 
vides a general, measurement-free scenario to control the 
h z noise by operating the NV as a magnetometer: any 
scheme in which the NV steady state P(ft) and hence 
the NV-induced steady-state nuclear polarization H(ft) is 
sensitive to the value ft, of h z could significantly suppress 
or amplify the h z noise. A possible scheme is to use very 
weak optical pumping at the NV ground state anticross¬ 
ing to polarize 13 C nuclei without significantly degrad¬ 
ing the NV sensitivity Aft ~ NV linewidth, ultimately 
limited to > I/T^nv by the true NV dephasing time 
± 2 ,nv ( n ot the inhomogeneous dephasing time T 2 NV ). 
For Ah <C ftmax, we can tune ft* to the region with the 
most negative H'(ft*) ~ — h max /Ah to reduce a from 
a th = ay/N /2 to an IV-independent value a ~ \J a Ah. 
For an estimate, we take ft max = 1 MHz to obtain the 
noise reduction factor = \/h max /Ah ss 10 (for 

Ah = 10 kHz) or 30 (for Ah = 1 kHz). 

The third approach to suppressing h z noise is (iii) to 
generate a sharp dip in the denominator of Eq. ©, e.g., 
by /^-dependent coherent population trapping (CPT) 
[43|. Now we show that this mechanism leads to the 
first observation of 13 C nuclear spin noise suppression 
0 - The setup of Ref. 0 consists of a A system (| ± 1) 
and lAi)) and a two-level system (|0) and | E y )), both un¬ 
der resonant pumping. Relevant processes are shown in 
Fig . [Ha) and the NV Hamiltonian can be found in Ref. 
[19( or [35[. The NV-nuclei coupling includes the contact 
HFI (AgSg + A e S e ) • io with the 14 N nucleus Io and the 
dipolar HFI J2n =l S • A„ • l n with the 13 C nuclei {!„}, 
where S g (S e ) is the NV ground (excited) state spin and 
S = S g + S e . The electron Knight fields on Io and I n are 
along e z || N-V axis and e z ■ A n = a n e n ,z, respectively. 
So we define the local coordinate (e n> x,e n ,Yi e n,z) and 
decompose the HFI into K = S StZ h z and the remain¬ 
ing part V, where h z = A g I 0iZ + a n /„ jZ . We define 
the nuclei basis |m) = |?n 0 ) ®n=i \ m n) as the product 
of eigenstates of each nucleus: Io, z \mo) = mo|mo) and 
I n ,z\n^n) = Wnjmn). Then we calculate the transition 
rates from Eq. (H|) (the 13 C nuclear spin flip by the 
NV ground state has no Fermi golden rule part, only 
the coherent part contributes) and solve Eq. Q 

numerically to obtain the steady-state nuclear spin pop¬ 
ulations {pm.m}. The intrinsic 13 C- 13 C interaction and 
13 C relaxation is included as a phenomenological depo¬ 
larization rate 7 c for each 13 C nucleus. The calculated 
14 N population on \mo = 0) and the relevant time scale 
~ 200 /xs (vs. experimental value ~ 353 ± 34 fj, s) agree 
reasonably with the experiment [Fig. QJb)]. We also con¬ 
firm that the observed decrease of the population at large 
H.A arises from the off-resonant excitation to IA 2 ), as ex¬ 
pected in Ref. 0 . An impressive observation [0 is 


the suppressed h z noise from 13 C nuclei in the absence 
of appreciable 13 C polarization, manifested as the nar¬ 
rowed CPT dip of the NV fluorescence. Using realistic 
and experimental parameters, we obtain the steady nu¬ 
clear populations {p m , m } at the preparation magnetic 
field and use them to calculate the unconditional and 
post-selected population on \E y ) at the readout magnetic 
field 0. When normalized to unity at large readout 
magnetic field B re = w re /(g e fiB), the results agree rea¬ 
sonably with the experimental fluorescence [Fig. [l](c)] . 

Finally we use Eq. (J9J to analyze qualitatively how 
the NV detects and suppresses h z noise in the CPT ex¬ 
periment [0. Given an instantaneous value h. of h z , the 
NV rapidly records h as a two-photon detuning in the 
NV steady-state P(h). Our calculation shows that the 
NV-induced nuclei flip always yields vanishing steady- 
state polarization M(h) = 0, so ft* = 0, a = <7 t h, and Eq. 
dU gives p(h) oc e - ft 2 /( 2 CT tii)[H r -|- < _j.(ft) +W(j. < _'f(ft)] -1 . The 
key is that the nuclei flip rates IU| < _ 4 .(ft), W\,^-\(h) oc \E y ) 
population (35 ], which has a sharp dip at the two-photon 
resonance h = 0. This generates a sharp peak in p(h) and 
hence suppresses the fluctuation of h z . Further analysis 
shows that at ft, = 0 , off-resonant excitation to IA 2 ) gives 
rise to non-vanishing \E y ) population and hence 13 C spin 
flip that fundamentally limit the noise suppression effi¬ 
ciency. To avoid this limitation, a possible scheme is to 
exploit the strain-induced non-vanishing (E y \S z \E y ) and 
hence the term oc \E y ){E y \(I++I-) of the dipolar HFI. In 
a magnetic field that quantizes 13 C nuclei along the N-V 
axis, a negative feedback is expected to significantly sup¬ 
press the h z noise without being limited by off-resonant 
excitation to |A 2 ). 


To summarize, we have presented a quantum theory 
for using an optically pumped NV center as a tunable 
bath to engineer various dissipative and collective nu¬ 
clear spin dynamics, as illustrated by the control of the 
nuclear spin dephasing, relaxation, and squeezing. It also 
reveals a general and efficient way to control the nuclear 
spin noise and clarifies the physical mechanism leading 
to the first observation of nuclear spin noise suppression 
[0 . Apart from NV centers, our theory can be readily 
applied to other quantum information platforms such as 
quantum dots and defect centers to engineer the nuclear 
spin dynamics for coherence protection | 0 . quantum 
feedback control [45| , and preparing non-classical nuclear 
spin states for quantum metrology 21 ] and information 


storage [23, [24]. It could also be extended to explore 


the crossover of the nuclei dissipation from Markovian 
to non-Markovian j4G], and from quantum Zeno to anti- 
Zeno effect [47] with gradually decreasing dissipation of 
the NV bath (by reducing the optical pumping strength). 
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Section I provides a detailed derivation of Eqs. (2)-(6) of the main text. Sec. II provides a perturbative, explicit expression 
for the nuclear spin transition rate VV p< m , which can be easily applied to a given experimental setup. Sec. Ill details the NV 
Hamiltonian, NV-induced nuclear spin transition rates, and calculation of NV fluorescence under the CPT condition [? ]. 


I. DERIVATION OF NUCLEAR SPIN EQUATIONS OF MOTION 


The starting point is 

,,, r „ . (Pm.nUX <Wfm.n) ., lrT -,„ , 

Pm,n(0 — -£m,nPm,n(0 — ! ^ — *(Ol|[V(f),p(f)]|n), 

. , ™ {Pm,n(0> ^^m.n) .„ , \ 

Pm,n(0 = -iTr e - 2 -iTr e <m|[V(0,p(f)]|n). 


(1) 

( 2 ) 


On the coarse grained time scale T e <$; At <s T\, T 2 , we treat „(fj) as slow variables, others as fast variables, and 6K mn and 
V(t) as first-order small quantities. Correspondingly, we decompose p(t ) = p (0) (0 + p (1) (0 + ■ ■ ■, where p (k \t) is a Ath-order small 
quantity. Since we always work in the nuclear spin interaction picture, the zeroth-order dynamics vanishes: (p)o = 0. 

The first-order dynamics (p m , n )i [Eq. (2) of the main text] is obtained from Eq. ([2} by replacing p(t) with p (l], (t) = 
Em.n l m )( n IPmji(0, obtained by coarse-graining 


Pm,n(0 — -£m,n/>m,n(0 


for an interval T e <sc At <sc 7], 72. Here P m n is the normalized electron steady state determined by -C m , n P m ,n = 0. The 
second-order dynamics (p m , n )2 is obtained from Eq. (0 by replacing p{t) with p (l \t), the solution to 


(Pm!n)l + Pm.'n = Em.nPm'n 


a(1) 


(1) . IPnijn hW m ,n! y „(())■, 

m'n - l -,-*<m|[V,p tU », 


where (pm?„)i = An,n(Pm,n)i comes from the first-order evolution (/? m n )i. Substituting into the above equation and coarse- 
graining for an interval \/y e <k Ar <sc 7), To gives 


p£ - 


) _ ■ r -i {Pm.n; d7f m n } . r \-Uf> f> Hr \ f> 

,11 — l P m.n-C m n " + l y \-L m.n " ' * in.p' pn ~ \ t nip/ pn I ni n f Pp,n 


f y 1 ^ -Tm.n + frUp n ) (Pm.pVp.n (Vp.n)m.pPm,n)Pm.p? 


where (X m , n + ioj mv ) 1 = - f 0 °° e (£m » +,Wm **dt, E,,, 1 ,, = - f Q °° e £m »'dt = lim^oC-Cm.n + M _1 , <• • • ) m .n = Tr e (- • • )P m ,n, and we 
have assumed V mn (t) — V mn e " JJ,nn ''. Then replacing pit) with p n> (t) in Eq. 0 gives the desired second-order nuclear spin 
dynamics. For m = n, neglecting the coupling of p m , m to other p p q (p + q), which amounts to neglecting the small second- 
order corrections to (Vitj) and electron-mediated nuclear spin interactions (they induce nuclear spin diffusion and depolarization, 
which can be included phenomenologically at the end of the derivation), we obtain Eq. (3) of the main text. The transition rate 


Wp«—m — 2 Re Tr e Vm,p(-C p ,m + hUp.m) (Vp.mPm.m \ I p.Ill /mm Pp.m )■ 


Using (-Cp.m + ia) p , m ) 1 P pm = Pp,m/(ito p ,m) and / J pm * P mm (|p) and |m) differs only by a single nuclear spin flip), we see that 
the term oc / J p m is (2/<y p , m ) Im(V m , p ) p . m (V p . m ) m , m « 0 and hence recover Eq. (4) of the main text. For m + n, neglecting the 
coupling of p mn to other variables and second-order self-energy corrections, which amounts to neglecting the small second-order 
corrections to (V(f)) and electron-mediated nuclear spin interactions, we obtain Eqs. (5) and (6) of the main text, where 


W p <_ m |n 


2 Re Tr f V m ,p(-C p ,n + ZWp.m) (V p . m P m .n \ 1 pm/milfpn)■ 
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II. PERTURBATIVE, EXPLICIT EXPRESSION FOR W fi m 

If <p|V(0|m) has multiple frequencies that are widely separated compared with the nuclei transition/dephasing rates, then 
I!p<—m is the sum of contributions from different frequency components. So we consider the case (p|V(f)|m) = Fe~ ,0>t with a 
single frequency. Equation (4) of the main text gives W pi - m = -2ReTr,,F i 6r/*7 J , where Q = (ho + F) 1 , P = P m m , and 
F = Tip.m. The key is to calculate Q perturbatively by dividing -£(•) = —i[H, •] + FA*) into the unperturbed part 

F\*) = -i[H d , r -- {|i> < , 'l> *1 + Z 7 “ m • m 


and the perturbation 


F n d ( .) = -i[H nd ,.] + 2 • 10, 

f±i 


where H d = (H nd ) is the diagonal (off-diagonal) part of H and T, = 2/7//- When ||_£ d + Ml » ||-C nd ||, an explicit 

expression for Wp^ m is obtained by replacing Q by Q d - Q d FA d Q d and using Q d \k)(j\ = {F d + MT* \k)(j\ = \k)(j\(i/zk,j ) with 
Zkj = 8^ - Ej - to - /(Tj. + T,- - 2jk,k6k,j)/2- F° r F — E fti Afi\f)(i\, we obtain W p< _ m as the sum of the golden rule part 


U'J 


V}j,Vf,iPi,i' 


Zfj 


(3) 


and the coherent part 


W coh = ? Tm 

v 'p<-m ^ 1111 


2 


t Tjnd y* 

n i'J' V W 

Zf,f 


nd \ 


v h H T, f 




VfjPw 




(4) 


where O, ,, s (i\d\i'). 


III. NV HAMILTONIAN, NUCLEAR SPIN TRANSITION RATES, AND NV FLUORESCENCE UNDER CPT 

A. NV Hamiltonian under CPT 


The NV Hamiltonian for the coherent population trapping (CPT) experiment has been discussed in [? ]. Here we reproduce 
it with greater detail using our own notations. The strain effect is neglected as we have confirmed that it produces a small 
influence on the nuclear spin polarization and noise suppression. The NV states of relevance include 3 ground triplet states |0) 
(energy 0), | ± 1> (energy D gs ), 6 excited triplet states \E y ), \E X ) (energy s Ex = s Ey ), \E X ), \E 2 ) (energy s El = s El ), \A X ) (energy 
fiAj), |Ai) (energy eaA, and one metastable singlet .S') (energy es). A linearly polarized laser with electric field V.\e~ u “' , / 2 + c.c. 
and frequency a>\ — e a - D gs resonantly excites the ground states | + 1) to the excited state |Ai) and, at the same time, off- 
resonantly excited | ± 1) to the excited state \AA) with detuning A - ea 2 - e a . Another linearly polarized laser with electric field 
E 2 e~ !W2, /2 + c.c. and frequency oh = s Ex resonantly excites |0) to \E y ). The relevant optical transition matrix element of the 
electric dipole operator d = -ex are (Ai|d| ± 1) = ±d a ^ E e ±/(2 V2), (A 2 |d| ± 1) = id a , E e ±/(2 V2), and (£ v |d|0) = d a . E ^ x t V2, where 
e± = e A ± ic y and d aE is the reduced matrix element of the electric dipole moment [? ]. Thus the laser coupling Hamiltonian 


HAD = 




& Al .-\ + ie uf ’&A 1 ,+i + ie " p o- Al -i)e 


,-UO\t 




■E w', 


e 2 o- Eyfi + h.c ., 


where we have defined <r, ; - s \i){j\, Q E = E 2 ■ <£’ v |d|0) = E 2x d a yJ V2, and FLac'^ = Ei • <A 1 |d| + 1) = d aE (E i, A + iEij)/(2 V2). 
Defining the bright state | b) = (e _, ^| + 1) - e" fl \ - 1))/ V2 and dark state | d) = (e~'^[ + 1) + e ,,s | - 1))/ V2, the laser coupling 
Hamiltonian simplifies to 


HAt) = %(d- Al , b + itTA 1 4)e~ ,u ' t + ^e-^'a E 
\r>. 2 




V2 


;,,o 


+ h.c.. 
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In the rotating frame \i//(t)) connected to the laboratory frame |i^iab(0) via |= e' Ho, \i/J\ a b(t)) with Hq = Se x (o'e x ,e x + o") + 
SeMeuEi + °~e 2 ,e 2 ) + + ( '' r A 1 ,A 2 ) + e s & s ,s + D gs (fr hl + o"-i,-i), the NV Hamiltonian is 


He — U> e S „ + Ao~A 2 ,A 2 


VT 




where we have included the ground state Zeeman term g e /dBBS z = u> e S z g of an external magnetic field B along the z (N-V) axis. 
We also include the dissipation of the NV center in the Lindblad form, including the spontaneous emission from |A ]) and IA 2 ) to 
|±1) (rate y/2), from \E y ) to |0) (rate y), and non-radiative intersystem crossing from |A 1 ) to |S) (rate y^), from IA 2 ) to |5) (rate 
y i2 ), from |S) to |0) (rate y^i), from lily) to | ± 1) (rate y ce ), and pure dephasing y^ for each excited state. 


B. NY-induced nuclear spin transition rate 


To calculate the NV-induced nuclear spin transition rates, we need to determine the NV steady state P(li) from -i[H e +S z h,P\ + 
£. e P(h ) = 0, where H e + S z g h - H e \oi e ^>s=u> e +h and 5 = cj e + li is the two-photon detuning. The calculation is straightforward and it 
turns out that relevant matrix elements of P(h ) can be expressed in terms of the \E y ) population (E y \P{h)\E v ), e.g., (S\P{li)\S) = 
(2yce/y s )(Ey\P(h)\Ey), (0\P(h)\0) = [1 + (y + 2 y ce )/W E ] < E y \P(h)\E y ), and (0\P(h)\E y ) = [i(y + 2y ce )/D.E](E y \P(h)\E y ), where 
We = Q 2 E /(2y ce + y + y^) is the transition rate between |0) and \E y ), Wa = El 2 A /(y + y s \ + y v ) is the transition rate between 
| ± 1) and |Ai), and 771 = y ce ly s \ . We treat the off-resonant | ± 1) —> IA 2 ) excitation perturbatively. Up to zeroth order, near the 
two-photon resonance 6 « 0 , we neglect 0(d A ) and high order terms and obtain 

(E v \P m (h)\E v ) » D 0 -^—r, 

S1+6 l 


where 


Do* 



1 


2 + 2m ^ + 1^’ 

r\ 2 11 Wa We 

mm _ Wa _ 

4 + + t At)' 


with 772 = y s /(y s + y ce ) and we have used 771 <s 1. We find that only \E y ), |0), and \d) are significantly populated since 
Ice ^ y s > y v 1 , y. Under saturated pumping We » y, the leading order correction from the off-resonant excitation to IA 2 ) is 


{E y \P (2 \h)\E y ) » -i-^-(1 - 2(E v \P <0 \h)\E v )), 

2r\\ y + y s i 

where Wa 2 - (0^/2)(y + y s 2 + y^)/A 2 is the off-resonant transition rate between | + 1) and IA 2 ) and we have used y s 2 <s y, yji- 
Now the | E y ) population (E y \P {0) (h)\E y ) + {E y \P^{h)\Ey) no longer vanish at 5 = 0. This limits the efficiency of l4 N nuclear spin 
cooling (in agreement with the expectation of Ref. [? ]) and 13 C nuclear spin noise suppression. 

Next we use Eqs. (0 and 0 [with / J m m replaced by P(li) and neglecting the dependence of ii |) m on p and m] to calculate 
the 14 N and 13 C nuclear spin transition rate. The transition rate W mo+ \^ mo {h) - W,„ 0 -i^ mo (h) for 14 N from |777o) to |7?7o ± 1) oc 
/q|/wo) is dominated by the following contributions from different NV transitions: A 2 x g (E y \P(h)\E y ) from |0) —» | ± 1) and 
A 2 e ('LfXf)( E y\P(h')\Ey) from \E y ) -> |/), where, / runs over A U A 2 ,E U E 2 states,** s (y + 2y ce ) / D 2 gs , and Xf = (1/4 )(T/ + 
Xip)/(.B e - £/) 2 , with Sf the energy of |/) in the laboratory frame. These transition rates differ from the phenomenlogical 
expression A e (E y \P\E y ) in Ref. [? ]. Similarly, the transition rate W mn± i<- mn of the 77 th 13 C nucleus from |m„> to \m„ ± 1) oc 
I„\m n ) is dominated by the following contributions: **(|A„__| 2 + |A„ j+ -\ 2 )(E y \P(h)\E y )/8 from |0) —> | + 1), lA^-l 2 ^, + 
XE t )(Ey\P(h)\E y )/2 from \E y ) -> |Ai),|£i), and \A„x.-\ 2 (xa 2 + XE 2 )(E y \P(li)\E y )/2 from \E y ) -> |A 2 ), \E 2 ), where A, hafi = e„, a • 
and e„ >± = &n,x — y • 


C. Calculation of NV fluorescence 

To compare with the experimentally observed suppressed 13 C nuclear spin noise, we first set a> e — 0.18MHz to obtain the 
steady nuclear state populations } and then calculate the nuclear-state-dependent NV population (E y \P mm \E y ) at the readout 
magnetic field w re = g,,fJi;B t - c . The fluorescence before post-selection is proportional to 2 m p m .m(E v \P m . m \E y ), while the fluores¬ 
cence after post-selection is proportional to Yum Pm,m{E y \P m ,m\E y )e~ eyTcmi ' E ^ Fmx '' [Ey \ where e is the photon collection efficiency 
and r cond is the photon collection time. 








